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1.  Introduction 


In  this  paper,  we  will  be  concerned  with  location  models  and,  in 
particular,  the  one  sample  location  model.  We  wish  to  make  inferences 
about  6,  the  center  of  symmetry  of  a  continuous,  symmetric  population 
with  c.d.f.  F,  based  on  a  random  sample  of  size  n.  Many  authors  have 
considered  this  problem,  although  many  place  more  restrictive  assumptions 
on  F,  such  as  normality.  If  these  additional  assumptions  are  questionable, 
practitioners  often  seek  procedures  which  remain  valid  when  no  specific 
form  for  the  underlying  distribution  is  assumed.  For  testing  hypotheses 
about  9,  the  Wllcoxon  and  sign  tests  have  emerged  as  two  of  the  more  popular 
procedures.  In  addition,  estimation  procedures,  both  point  and  interval, 
based  on  these  rank  statistics  have  been  developed.  However,  the  sign  test, 
while  optimal  if  F  is  the  double  exponential  c.d.f.,  is  quite  inefficient 
at  the  normal  model.  Conversely,  the  Wilcoxon  test,  which  is  optimal  if 
F  Is  the  logistic  c.d.f.,  is  less  efficient  for  heavier  tailed  distributions. 
Also,  the  estimation  techniques  derived  from  the  Wilcoxon  statistic  (Lehmann 
(1975),  Chapter  4,  Section  3)  are  computationally  difficult  for  moderate 
or  large  samples  (Johnson  and  Ryan  (1978)). 

Thus,  we  propose  a  class  of  rank  statistics  "between"  the  sign  and 
Wilcoxon  statistics.  The  corresponding  tests  and  estimation  procedures 
will  retain  much  of  the  simplicity  of  those  based  on  the  sign  statistic 
but  with  improved  efficiency  (especially  near  the  normal  model).  In  addi¬ 
tion,  all  the  procedures  have  excellent  robustness  properties. 


r 


2. 


2.  Testing 


Let  X^  £  ...  £  Xjj  be  the  ordered  values  of  a  random  sample  from  F. 

Let  0  <  a(l)  <  a(2)  <  ...  <  a(N)  be  a  given  set  of  scores.  Then  W„(0  )  ■ 

Za(i)V  (0  )  is  a  linear  signed  rank  statistic  for  testing  H  :  0*0  where 

1  o  o  o 

V.  (0  )  «  1  if  Xj  -  0  >0  and  0  otherwise  and  |X,  -  0  I  <  ...  <  [X,  -  0  I 

l  o  o  1  d^  o  —  —  1  d^  o 

In  this  case,  i  is  the  rank  of  |X^  -  0q |  among  the  absolute  values,  i  *  1, 

...,  n.  Notice  that  a(i)  *  i  yields  the  Wilcoxon  statistic  and  a(i)  *  1 

yields  the  sign  statistic.  Our  approach  is  to  select  a  set  of  scores 

"between"  the  sign  and  Wilcoxon  scores.  Let 


C  k  *■  “  bk  +  •  •  * »  N 

where  0  «  yq  £  yL  <  y£  <  ...  <  yfc  £  yw  -  1,  bt  -  [Nyt]  t  -  0,  . . . ,  k  +  1 

and  [•]  is  the  greatest  integer  function.  Let  TN  k(0)  *  Zak(i)Vi(0).  Then 

T„  ,  is  called  the  k-step  statistic. 

N,k  r 

Noether  (1973)  studied  the  1-step  statistic  with  particular  emphasis 
on  interval  estimation  based  on  T^  Notice  that  if  y^  *  1  t*ie 

sign  statistic.  Policello  and  Hettmansperger  (1976)  discuss  the  1-step 
statistic  with  regard  to  adaptive  methods  of  estimation  of  0.  Notice  that 
the  freedom  to  vary  y^,  ...,  y^  makes  this  a  rich  family  of  statistics  and 
hence  some  member  of  the  family  should  be  suitable  for  testing  Hq  over  a 
wide  range  of  symmetric  underlying  distributions. 


3. 


To  test  H  :  0*8  vs.  H  :  0  >  0  ,  we  will  reject  H  if  T„  .  (8  )  >  c. 

o  o  a  o  o  N,k  o 

However,  we  will  consider  only  the  case  8  *0  since  otherwise  we  can 

o 

subtract  0  from  each  sample  observation.  Henceforth,  we  denote  TM  ,  (0) 

®  Nf  K 

by  T„  ,  .  Next,  we  consider  the  null  distribution  of  T„  ,  . 

*  JC  N,  K 

k 

Theorem  1.  If  0  «  0,  then  TXI  .  ■  £  iW.  where  W.  is  distributed  as  binomial 

NjK  ^ 

with  parameters  bi+1  -  bt  and  1/2,  and  . Wfc  are  independent. 

This  result  follows  directly  from  Lemma  10.1.11  of  Randles  and  Wolfe 
(1979). 

Corrollary.  If  8*0,  then 

E<Yk>  •  X  -  b!>/2-  T"<Yk>  ■/,  l2<bi+i  -  V'4 

i*l  i*l 

Notice  that  T^  ^  follows  a  binomial  distribution  with  parameters  N  -  b^ 

and  1/2  and  that  T  _  is  a  weighted  combination  of  two  independent  binomial 

random  variables.  We  conclude  discussion  of  the  null  distribution  of  T„  , 

N,k 

with  the  following  result  on  the  asymptotic  distribution  of  T  ,  . 

y  K 

Theorem  2.  Under  Hq:  0*0,  for  k  >_  1  fixed,  the  limiting  distribution 

(as  N  -*■  «)  of  (TN>k  -  E(TNjk))/(Var(TN>k))1/2  is  N(0,  1). 

Proof.  Theorem  1  exhibits  T„  ,  as  a  linear  combination  of  a  fixed  number 
-  N,k 

of  independent  binomial  random  variables.  Since  each  of  these  converges 
to  a  normal  distribution,  the  result  follows. 

Example.  Thirty  observations  were  generated  from  a  normal  distribution 
with  u  *  0,  a  *  4.  The  ordered  data  (rounded  to  two  decimal  places)  follows. 


-10.33 

-2.83 

0.01 

2.30 

4.28 

-6.16 

-1.36 

0.05 

2.38 

4.30 

-6.03 

-1.30 

0.31 

2.97 

4.92 

-5.80 

-0.93 

1.25 

3.47 

6.19 

-5.10 

-0.53 

1.81 

3.67 

6.32 

-3.71 

-0.52 

2.24 

4.08 

7.96 

4. 


We  will  compute  the  2-step  statistic  with  steps  at  y^  *  .2,  y^  *  .6.  The 

reason  for  this  choice  will  be  apparent  later  In  this  section.  Since  N  - 

30,  *  [30(.2)]  -  6,  bj  ■  18.  Thus,  ^  Ignores  the  six  observations 

with  smallest  absolute  value,  counts  the  number  of  positive  observations 

with  absolute  values  ranking  from  7  to  18,  and  adds  to  this  sum  twice  the 

number  of  positive  observations  with  absolute  values  ranking  above  18. 

Thus,  T„  .  *  8  +  2(7)  *  22.  For  testing  H:  8*0  vs. H:  0  >  0,  we  have, 

« o  a 

1/2 

correcting  for  continuity,  Z  ■  (T„  .  -  .5  -  E(T„  -))/(VarTu  _)  and  so 

n9k  n,z  n , z 

Z  -  (22  -  .5  -  18)//L5  -  .90  or  p  -  .1841. 

The  exact  distribution  of  T^  under  Ha  is  quite  complex  and  will  not 

be  discussed.  We  mention  only  that  Markowski  (1980)  has  shown  that  the 

exact  distribution  of  T„  .  depends  only  on  P(X.  +  X4J,  >  0)  j  -  1,  .... 

n,j.  j  j+oi 

N  -  b^,  but  the  form  of  the  dependence  is  quite  complex. 

In  order  to  discuss  asymptotic  properties  of  the  inference  procedures 
studied,  it  is  necessary  to  discuss  some  aspects  of  the  asymptotic  distri¬ 
bution  of  TN  k  under  Ha>  First,  Markowski  (1980)  has  shown  that  TN  ^  is 
asymptotically  normally  distributed  under  a  sequence  of  contiguous  location 
alternatives.  The  proof  involves  a  modification  of  a  technique  employed 
by  Albers,  Bickel,  and  van  Zwet  (1976).  Next,  when  the  limit  exists,  define 

the  Pitman  efficacy  of  a  statistic  Rjj  by  C(IO  =  lim  E©^)!  /  (N  VarQ  R^)1^ 

0*0 

where  the  variance  is  computed  under  the  null  hypothesis  and  the  mean  is 
computed  under  the  alternative.  Under  regularity  conditions,  the  asymptotic 
performance  of  two  test  statistics  can  be  compared  by  considering  the  ratio 
of  their  squared  efficacies.  See  Lehmann  (1975,  Appendix,  Section  6)  for 
a  complete  discussion. 

The  efficacy  of  the  k-step  statistic  can  be  computed  directly  from 


the  following. 


Theorem  3.  If  F  has  finite  Fisher's  information,  and  the  density,  f,  of 
F  is  differentiable,  then 


c2(IN,k> 


4 [/  J+(2F(x)  -  1)  f'(x)dx]2 

o 

I 

f  J+(u)du 
o 


(2.2) 


where  J+,  the  limiting  scores  function,  is  given  by 


0  0  £  u  <  y 

1  Yx  <  u  <  Y2 


J+(u,  k) 


k  Yk  <  u  <  1 


(2.3) 


Proof.  Begin  with  10.2.11,  p.  338  of  Randles  and  Wolfe  (1979)  and  apply 
an  integration  by  parts  to  yield  the  result. 

Corollary.  Under  the  conditions  of  Theorem  3,  the  squared  efficacies  of 
the  one  and  two  step  statistics  are: 


C2(T..  .)  .  +  rp/2]) 

1  -  Y, 


(2.4) 


and 


c2<V> 


[f{F_1[(l  +  y1)/2]}  +  f{F-1((l  +  Y2)/2]}]2 
1  -  Y,  +  (Y2  -  Yx)/4 


(2.5) 


Using  (2.4),  we  determine  in  Table  1  that  member  of  the  class  of  one- 
step  statistics  which  maximizes  the  efficacy  for  the  following  underlying 
distributions: 

A1:  Normal,  f[F_1(x)]  -  (2TT)"1/2{exp{ (F_1(x))2/2]}, 

A2:  Logistic,  f[F-1(x)]  -  x(l  -  x) , 


A^:  Contaminated  Normal,  e  •  .05,  a  *  9,  f(x)  -  (1  -  e)<J>(x)  +  ^(x/a) 
where  $(*)  is  the  standard  normal  density  function, 

A^ :  Double  Exponential,  f[F-1(x)]  *  x  if  0  <  x  <  .5  and  1  -  x  if  .5  <  x  <  1. 

-  Table  1  about  here  - 

Table  1  indicates  that  for  the  1-step  statistic,  the  optimal  choice  of 
y^  decreases  as  the  tails  of  the  underlying  model  become  heavier.  Also,  we 
see  that  despite  the  simplicity  of  the  1-step  statistics,  the  best  1-step 
statistic  is  quite  efficient  compared  to  the  optimal  choice  in  the  class  of 
possible  statistics.  In  addition,  we  have  included  a  comparison  of  the  best 
1-step  statistic  as  well  as  the  1-step  with  y  *  1/3  as  opposed  to  the  sign 
statistic.  These  latter  two  comparisons  were  first  discussed  by  Noether 
(1973)  in  a  somewhat  different  setting.  Notice  that  the  statistic  with 
*  1/3  seems  to  be  an  improvement  over  the  sign  statistic  since  its 
efficacy  is  more  stable  and  improves  significantly  over  the  sign  statistic 
at  the  normal  model. 

-  Table  2  about  here  - 

Table  2  illustrates  the  effect  on  C^(T„  „)  caused  by  changes  in  y, 

N ,  l  'X 

and  Y2*  addition,  the  optimal  choice  of  (y^  Y2)  is  noted  for  each 
underlying  model.  Notice  that  the  efficacy  changes  only  slightly  near  the 
optimum  value  so  that  a  choice  of  (y^,  y^)  in  a  neighborhood  of  the  optimal 
value  will  still  yield  nearly  highest  efficiency.  Of  course,  typically  the 
underlying  distribution  may  not  be  known  and  so  a  test  with  good  efficiency 
over  all  distributions  considered  would  be  desirable.  From  Table  2,  we  see 
that  the  choice  y^  ■  .2,  Y2  “  *6  has  such  properties.  In  Table  3,  we  compare 
this  2-step  statistic  with  the  sign,  Wllcoxon,  and  t  statistics  as  well  as 


7. 


Che  likelihood  ratio  statistic  for  the  particular  underlying  distribution. 

We  see  that  ^  improves  upon  the  sign  statistic  while  retaining  much  of 
its  simplicity.  In  addition,  it  compares  quite  favorably  with  the  t- 
statistic  at  the  normal  model  (e  -  .89)  while  being  superior  to  the  t-statistic 
at  the  other  distributions  considered.  We  see  also  that  this  2-step  statistic 
has  properties  very  similar  to  that  of  the  Wllcoxon  statistic. 

-  Table  3  about  here  - 

Thus,  if  one  step  is  used,  we  recommend  y^  *  1/3  and  if  two  steps  are 
used  take  *  .2,  y^  *  .6.  These  recommendations  are  based  on  the  effi¬ 
ciency  calculations  of  Table  3.  In  the  remainder  of  the  paper,  we  will 
investigate  the  estimation  procedures  associated  with  these  tests  and  study 
their  robustness  properties. 


8. 


Point  Estimation 


Hodges  and  Lehmann  (1963)  illustrate  a  method  for  deriving  point 

estimators  from  rank  statistics  such  as  We  denote  by  ek(Yr  ....  Yk) 

% 

or  simply  by  0  the  Hodges  Lehmann  estimator  (H-L)  of  0  derived  from  ^ . 

Bauer  (1972)  has  shown  that  the  linear  rank  statistics  W„(9)  ■  Ea(i)V.(9)  as 
defined  in  Section  2  are  non- increasing  step  functions  which  can  have  steps 
only  at  the  Walsh  averages,  {(x.^  +  Xj)/2  1  1  £  1  £  j  £  N).  The  sign  statistic 
has  steps  only  at  (x^  :  1  £  i  £  N)  while  the  Wilcoxon  has  steps  at  each  of 
the  Walsh  averages.  Following  the  method  of  Bauer  (1972),  we  find  that 
has  steps  of  size  1  at  the  averages  in  =  { (xi  +  x^)/2  :  j  -  i  =  bt»  t  *  1, 

.  ,  k  define  D. 


...,  k)  and  hence  0K  *  med  D^.  Now  for  t  =  1, 


k,t 


{(x.^  +  xb  +i)/2  :  i  ■  1,  ...,  N  -  bt>,  the  Walsh  averages  corresponding  to 

the  t'th  step  of  T  .  .  Then  D,  »  U  D  and  so  0k  *  med  (U  D,  ).  The 

advantage  of  this  representation  is  that  the  Walsh  averages  in  ^  are 

a  priori  ordered.  This  means  that  for  small  k,  finding  the  median  0  may 

not  be  difficult.  For  example, 

f  y  +  x  }  /2  if  N  -  b.  odd 

iX(N  -  b,  +  l)/2  X(N  +  b,  +  l)/2'  1 

02  -  1  (3.1) 

{X(N  -  b1)/2+  X(N  +  b1)/2  +  X(N  -  \  +  2)/2  + 

X(N  -  b2  +  2)/2^/4  if  N  "  bi  even* 

Thus,  9  can  be  computed  directly  from  the  ordered  sample.  As  another 
example,  we  use  the  ordered  structure  of  the  sets  and  to  exbibit 
a  simple  algorithm  for  computing  9  .  Our  approach  is  to  begin  with  medians 
of  each  of  the  sets  and  and  then  by  a  series  of  comparisons  of 

A  A 

Walsh  averages,  determine  9  . 
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Let  “ij  *  (xi  +  Xj)/2  ***  B(a’  b)  “  ma,  b1+a  -  “b,  b2+b 

a  *  1 . N  -  b^  b  *  1,  N  -  bj 

Algorithm 


Case  1.  N  -  b^  even,  N  -  b2  odd. 

First  consider  B((N  -  b^)/2,  (N  -  b2  +  1)1 2)  and  suppose  for  now  that 

it  is  greater  than  0.  Notice  that  this  implies  m,„  ,  w_  <  0  < 

v,N",D2^/^»  V.N+D2; / 2 

“(N-b^/^,  (N+b^/2.  Next,  compute  B((N  -  b.^/2,  (N  -  b2  +  l)/2  +  1).  If 
B((N  -  bx)/2,  (N  -  b2  +  l)/2  +  1)  <  0,  then  0  =  ^.h^)/2t  (N+b^/2'  If 

B((N  -  b1)/2,  (N  -  b2  +  l)/2  +  1)  >0,  compute  B((N  -  b^/2  -  1,  (N  -  b2  +  l)/2 


+  1). 


If  B((N  -  bx)/2  -  1,  (N  -  b2  +  l)/2  +  1)  >  0,  compute 


B((N  -  b^)/2  -  2,  (N  -  b2  +  l)/2  +  2).  Continue  until  B((N  -  b^)/2  -  i, 

(N  -  b2  +  1 ) / 2  +  i)  >  0  and  B((N  -  b^/2  -  (i  +  1),  (N  -  b2  +  l)/2  +  (i  +  1)) 
<  0.  Then  if  B((N  -  b^/2  -  i,  (N  -  b2  +  l)/2  +  (i  +  1))<  0,  9  = 

m(N-b1)/2-i,  (N+h1)/2-i  and  if  B((N  *  bl)/2  ~  i’  (N  “  b2  +  1)/2  +  (i  +  1) 


>  0,  9  =  m 


(N-b 2+l ) / 2+ (i+1 ) ,  (N+b 2+l ) / 2+(i+l ) ‘ 


Recall  that  we  assumed 


B((N  -  b^)/2,  (N  -  b2  +  l)/2)  >  0.  If  not,  then  proceed  in  a  completely 


analogous  fashion  except  increase  the  first  component  and  decrease  the  second 


until  B((N  -  b^/2  +  i,  (N  -  b2  +  l)/2  -  i)  >  0. 

Case  2.  N  -  b^  odd,  N  -  b2  even. 

Proceed  as  in  Case  1  but  begin  with  B((N  -  b^  +  l)/2,  (N  -  b2)/2). 

Case  3.  N  -  b^  odd,  N  -  b2  odd. 

Begin  with  B((N  -  b^  +  l)/2,  (N  -  b2  +  l)/2).  Proceed  as  in  Case  1  reaching 
the  point  where  (e.g.)  B((N  -  b^  +  l)/2  +  i,  (N  -  b2  +  l)/2  -  i)  <  0  and 
B((N  -  bx  +  l)/2  +  (i  +  1),  (N  +  b2  +  l}/2  -  (i  +  1))  >  0.  If 
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B((N  -  b1  +  l)/2  +  (i  +  1),  (N  -  b2  +  l)/2  -  i)  >  0  (3.2) 

B((N  -  b1  +  l)/2  +  i,  (N  -  b2  +  l)/2  -  (i  +  1)  <  0,  (3.3) 

A 

then  0  -  Cm (N-b ^+1 ) / 2+  ( i+1 ) ,  (N+b^D/Z+U+l)  +  m(N-b2+l)/2-(i+l)  ,  (N+b2+l)/2- 

(i+l)}/2.  If  the  signs  of  the  inequalities  (3.2)  and  (3.3)  are  other  than 

/\ 

those  listed,  then  once  again  0  can  be  determined  in  an  analogous  fashion. 

Case  4.  N  -  b^  even,  N  -  b2  even. 

Begin  with  B((N  -  b^)/2,  (N  -  b2)/2)  and  proceed  as  in  Case  1.  Then  make 
comparisons  as  in  (3.2)  and  (3.3).  0  is  determined  by  the  direction  of 
these  inequalities.  Notice  that  it  may  happen  that  the  final  stage  cannot 
be  reached  because  the  second  component  can  no  longer  be  increased  (or 

/x 

decreased).  At  this  point,  0  can  be  computed  immediately  since  in  this 
case: 

1)  0  must  come  from  the  set  D^. 

2)  The  position  of  the  largest  element  of  D^2  in  the  ordering  of 

will  have  been  determined  by  previous  steps. 

Next,  we  illustrate  the  use  of  the  algorithm  using  the  data  of  the 

/x 

example  in  Section  2.  We  will  compute  0(.2,.6)  for  this  data.  Now  N  =  30, 
b^  *  6,  b2  =  18  and  so  N  -  b^  *■  24,  N  -  b2  =  12.  Thus,  we  seek  the  18th 
and  19th  largest  elements  of  .  To  find  these  we  apply  Case  4  of  the 
Algorithm.  We  begin  with 


B(12,6) 


X12+X18  X6+X24 

2  "  2 


1.35  >  0. 


Next  we  compute 


B(ll, 7) 


*11*17 

2 


X7+X25 


-0.17  <  0. 


Thus,  we  have  reached  the  stopping  point  and  @(.2,. 6)  can  be  determined 
after  two  more  comparisons. 


Now  B(12, 7)  -  .27  >  0,  B(12,8)  -  -1.22  <  0,  B(ll,6) 


.91  >  0.  Hence,  the 


7th  smallest  element  of  is  between  the  llc^  and  12 
and  so  It  Is  18tl*  smallest  in  D^.  Similarly,  the  12C^ 


D31  is  19 


th 


smallest  in  D^.  Thus 

X  +X  X  +X  g 

0(.2, .  6)  +  2_^.  ] 


.79 


smallest  in 
smallest  element  of 


Note  that  for  this  data  X  «  .46,  the  sample  median  is  .78,  and  the  median 
of  the  full  set  of  Walsh  averages  is  .705.  The  mean  seems  to  have  been  greatly 
influenced  by  the  observation  -10.33.  If  we  remove  this  point,  the  sample 
mean  becomes  .837,  which  is  similar  to  the  other  estimators.  Notice  also 
that  5  Walsh  averages  and  5  comparisons  were  needed  in  the  computations. 

This  can  be  contrasted  with  the  465  Walsh  averages  which  enter  into  compu¬ 
tation  of  the  median  of  the  full  set  of  Walsh  averages,  which  is  the  estimator 
based  on  the  Wllcoxon  statistic. 

For  k  >  3,  the  same  technique  can  be  applied  although  the  number  of 
ordered  sets  increases  and  as  a  result  substantially  more  averages  may  be 
needed  before  the  estimator  can  be  determined. 

A  A 

Next,  we  consider  the  efficiency  of  0  .  Let  p^  and  be  two  consistent 

estimators  of  6.  Suppose  that  N172  (p^  -  0)  has  a  N(0,a2(Pi>)  limiting 

distribution.  Then  the  asymptotic  relative  efficiency  of  P^  with  respect 

to  p^  is  ARE(p1,p2)  *  a2(p2)/a2  (p^) .  In  addition,  if  p.^  is  the  H-L  estimator 

2  ~  2  -1 

derived  from  a  statistic  R^,  then  a  (p^)  =  [C  (Rjj)]  and  hence  the  ARE  of 
two  H-L  estimators  is  the  same  as  the  Pitman  efficiency  of  the  corresponding 
test  statistics.  See  Randles  and  Wolfe  (1979,  p.  227).  In  particular, 
since  TN  k  is  asymptotically  normal,  JT172#  -  0)  has  N(0,[C2(TN>k)]_1)  as 
its  limiting  distribution,  where  IC2(TN  k>]  1  is  given  in  (2.2). 


A 
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As  a  result,  the  efficiency  comparisons  of  Table  1-3  continue  to  hold 

when  comparing  the  corresponding  point  estimators.  In  particular,  using 

Table  3,  we  can  compare  6  (.2,. 6)  with  the  sample  mean,  median,  and  the  median 

of  the  set  of  all  Walsh  averages.  We  see  that  6  (.2,. 6)  compares  favorably 

with  each  of  these  estimators,  and,  in  particular,  its  efficiency  very  closely 

parallels  that  of  the  median  of  all  Walsh  averages.  However,  computation  of 
~3 

6  (.2,. 6)  requires  only  a  small  subset  of  all  Walsh  averages,  as  we  have 

-'3 

previously  seen.  Thus,  0  (.2,. 6)  is  a  viable  alternative  to  the  other  esti¬ 
mators  if  both  simplicity  and  efficiency  are  considered. 

Next,  we  propose  an  alternative  estimator  of  0  which  further  utilizes 
the  ordering  of  the  elements  of  and  D^*  The  estimator  suggested  is  a 
linear  combination  of  the  medians  of  the  sets  and  where  the  weights, 

A,  and  1  -  A,  are  proportional  to  the  number  of  elements  in  and  D^. 
respectively.  Following  this  approach,  we  obtain 

9*(Yl»Y2)  *  (X[(l-y  )N/2]  +  X[(l+y  )N/2])A/2  + 

1  (3.4) 

(X[(l-Y2)N/2]  +  XI(1+y2)N/2])(1  '  X)/2 

A 

where  A  ■  proportion  of  elements  in  which  are  in  D^.  Although  0*  is  not 

a  H-L  estimator  based  on  any  member  of  T„  _ ,  we  would  expect  it  to  have 

N,2 

A 

similar  properties.  This  is  confirmed  in  Table  4  which  compares  0*(.2,.6)  = 

s'  ’ 

^X[.4N]  +  X[  6N]^3  +  ^X[.2N]  +  X[.8N]^6  W±th  e3(’2’*6)  as  wel1  as  with 
other  commonly  used  estimators  of  0.  We  see  that  very  little  efficiency  loss 

A  /V 

results  from  using  0*  in  place  of  0,  and  some  computational  advantage  is 
gained.  Again  note  the  improved  efficiency  of  this  simple  estimator  at  the 
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normal  distribution.  Thus,  in  terms  of  simplicity  and  efficiency,  both 

A 

9*(.2, .6)  and  9(.2,.6)  are  viable  alternatives  to  the  classical  estimators 
considered  in  this  study. 

-  Table  4  about  here  - 

We  conclude  this  section  vith  a  brief  discussion  of  the  favorable 

A  Al. 

robustness  properties  of  9  and  9*.  The  influence  curve  of  9  can  be 
determined  using  the  methods  of  Hettmansperger  and  Utts  (1977)  and  that 
of  9*  by  following  Huber  (1977).  In  each  case,  we  find  that  the  Influence 
curve  is  a  step  function  with  a  finite  number  of  jumps  so  that  it  is  bounded 
but  discontinuous.  The  boundedness  is  a  desirable  property  since  it  indi¬ 
cates  that  the  estimator  is  not  unduly  influenced  by  any  single  observation, 
regardless  of  its  magnitude.  The  discontinuity  of  the  influence  indicates 
instability  in  the  estimator  at  the  points  of  discontinuity.  A  similar 
effect  has  been  noted  by  Hampel  (1974)  for  the  Winsorized  mean.  This  in- 

A 

stability  in  9  is  further  illustrated  for  the  case  k  =  3  by  the  following 

~3  — 

result  which  gives  a  worst  case  comparison  of  9  with  X. 

Theorem  4.  Let  ft  be  the  class  of  symmetric,  continuous,  unimodal  densities. 
Then 

inf  ARE(9(y  ,y2),  X)  *»  (1  -  Y1)3/(3d  -  Y2)  +  (Y2  -  Y ^M)  (3.5) 

feft  '*  1 

We  defer  the  proof  until  the  appendix.  Notice  that  when  y^  *  *2  and  y2  ■ 

.6,  the  infimum  is  .39.  This  is  a  slight  improvement  over  the  median  vs. 

the  mean  for  which  the  infimum  is  .33.  The  infimum  in  the  case  of  the 

median  of  the  Walsh  averages  vs.  the  mean  is  .864.  In  addition,  it  should 

be  noted  that  the  least  favorable  distribution  is  one  which  exhibits  maximum 

Instability  at  those  points  which  are  discontinuities  of  the  influence  curve. 


4« 
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Another  measure  of  robustness  Is  the  breakdown  point  (Huber  1981) . 

While  the  influence  curve  measures  the  effect  of  an  additional  observation 
at  y  on  the  estimator,  the  breakdown  point  gives  the  fraction  of  gross 
errors  the  estimator  can  tolerate  before  it  becomes  unbounded  -  Hampel 
(1971)  formulates  a  technical  definition  as  well  as  some  examples.  The 
breakdown  point  of  0*(y1>  y2)  is  (1  -  Y2)/2.  This  follows  from  Huber  (1981) 
but  is  also  intuitively  clear  since  (1  -  y2)/2  is  the  most  extreme  quantile 

A 

used  to  compute  9*.  Following  Hettmansperger  and  Utts  (1977),  it  can  be 
shown  that  the  breakdown  of  0  (y^,  y^)  is 


b(93(Y1,  y2))  - 


(Y2  -  Y1)/2 
2"1  -  (Yx  +  Y2)/4 


3Y2  1  Yx  +  2 
3y2  <  y1  +  2 


Notice  that  the  breakdown  of  0*(.2,.6)  is  .2  while  that  of  0(.2,.6)  is  .3. 

_  <v 

We  also  note  that  the  breakdowns  of  X,  HL,  and  X  are  0,  .293,  and  .50, 
respectively.  Table  5  indicates  how  the  breakdown  of  Q(y^t  Y2)  varies  with 
Y-L  and  y2- 


-  Table  5  about  here  - 


4.  Interval  Estimation 


Confidence  Intervals  for  6  based  on  rank  tests  are  described  in  Lehmann 

(1963),  Sen  (1966),  and  Bauer  (1972).  Applying  the  methods  of  Bauer  to  the 

A|/ 

statistics  we  see  that  the  endpoints,  0^  and  9y,  of  the  confidence 

interval  must  be  elements  of  D^.  In  particular,  if  p0.q(tn  k  <_  C)  »  a/2, 

A|f 

then  a  (1  -  a)  100%  confidence  interval  for  9  has  endpoints  9  ,  the  (C  +  l)'st 

L 

ordered  element  in  and  9^,  the  (N^  -  C) 'th  ordered  element  in  Dfc,  where 

N,  ■  kN  -  Z  b^,  the  number  of  elements  in  D,  .  Thus,  the  problem  of  computing 
i-1  * 

the  endpoints  of  this  interval  is  reduced  to  that  of  finding  particular 
Walsh  averages  within  the  restricted  set  D^. 

If  k  ■  2,  we  have  previously  seen  that  the  elements  of  D2  are  ordered 
and  so  we  have  0L  *  (xc+1  +  Xc+l+b1^2’  and  ®u  “  ^N-bj-C  +  ^-C^2  where 
C  can  be  determined  using  the  fact  that  2(0)  'V  Bin(N  -  b^  1/2)  or  by  a 
normal  approximation.  This  result  was  first  obtained  by  Noether  (1973). 

If  k  *  3,  we  seek  the  (C  +  l)'st  and  (2N  -  b^  -  b2  -  C) 'th  ordered  Walsh 
averages  in  D^.  The  algorithm  of  Section  3  modified  to  start  from  selected 
quantiles  of  the  sets  and  can  be  used  to  compute  these  endpoints. 

We  illustrate  this  with  the  data  previously  considered.  We  will  calculate 
an  approximate  95.44%  confidence  interval  for  6.  Using  the  Corollary  to 
Theorem  1,  a  normal  approximation  with  continuity  correction  yields 
C  i  <T1(a/2)((b2  -  bj_)/4  +  N  -  b2)1/2  +  (b2  -  bx)/2  + 


N  -  b2  -  1/2. 


(4.1) 


Now  letting  “  -2,  y2  *  -6  and  recalling  that  N  «  30,  we  obtain  C  *  10. 

A  A 

Thus,  0L  is  the  11 'th  smallest  element  in  and  9y  is  the  26 'th  smallest 
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element  In  Dy  To  compute  9^,  notice  that  the  number  of  Walsh  averages  In 
0^  is  12.  As  a  result,  we  first  consider  B(7,4)(since  ll(24)/36  ■  7  and 
ll(12)/36  A  4).  in  applying  this  method,  care  is  needed  to  assure  that  the 
starting  point  be  of  the  form  B(a,b)  with  a  +  b  ■  d  if  the  d'th  smallest 
element  is  desired.  Now,  B(7,4)  <  0  so  we  check  B(8,3)  >  0.  Furthermore, 
B(8,4)  >  0  and  B(7,3)  >  0.  Thus,  the  7'th  element  of  is  ll'th  smallest, 
i.e.  0^  »  (X^  +  X^)/ 2  ■  -1.41.  By  a  similar  argument  starting  with  B(17,9) 
we  find  B(17,9)  >  0,  B(16,10)  <  0,  B(17,10)  >  0,  and  B(16,9)  >  0.  Thus, 

§U  *  <*16  +  X22  )/2  «  2.36.  Notice  that  only  8  Walsh  averages  were  required 
in  this  case,  whereas  there  are  a  total  of  465  Walsh  averages  for  this 
sample.  The  approximate  95.44%  confidence  interval  based  on  the  Wilcoxon 
statistic  which  requires  ordering  of  these  465  Walsh  averages  is  (-1.04, 

2.30)  for  this  data. 

Lehmann  (1963)  has  considered  two  definitions  of  the  asymptotic  relative 
efficiency  of  two  interval  estimation  procedures;  one  in  terms  of  the 
probability  of  covering  a  false  value,  the  other  in  terms  of  the  asymptotic 
behavior  of  the  lengths.  Randles  and  Wolfe  (1979)  denote  these  by  ARE  and 
L-ARE,  respectively.  Lehmann  has  shown  that  in  many  regular  cases,  including 
all  those  considered  in  this  work,  the  two  definitions  yield  the  same  effi¬ 
ciencies.  Furthermore,  the  efficiencies  of  confidence  intervals  derived 
from  test  statistics  using  the  method  of  Lehmann  is  just  the  ratio  of  squared 

•'vJC 

efficacies.  Hence,  we  see  that  the  favorable  efficiencies  of  T^  and  9 
also  apply  to  the  interval  (0^,0^).  In  particular,  if  Y^  ■  .2,  Yj  ■  .6, 

*  3  /v  3 

the  interval  (0T  ,0  )  has  efficiency  similar  to  the  interval  derived  from 

la  U 

the  Wilcoxon  statistic,  yet  requires  computation  of  only  a  very  small  number 
of  Walsh  averages. 
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Next,  we  consider  another  confidence  interval  for  9  obtained  by  deter¬ 
mining  two  separate  intervals,  one  from  and  the  other  from  D32*  and  then 
combining  the  endpoints  of  these  intervals.  Suppose  thaC  the  desired  con¬ 
fidence  coefficient  is  1  -  a.  Then  an  approximate  (1  -  a)  1002  confidence 

interval  for  0  based  on  the  elements  of  D,,  will  be  (W_  , . ,  W„  _  )  where 

Jl  C^+l  N-b^-C^ 

W,  <  W.  <  ...  <  W„  ,  are  the  ordered  elements  of  D,.  and  using  a  normal 
1  —  1  —  —  N-b.  Ji 

—  1  1  /2 

approximation  $  (a/2)  (N  -  b^)  /2  +  (N  -  b^  -1)  /2.  Similarly, 

an  approximate  (1  -  a)  100%  confidence  interval  for  9  based  on  the  elements 
of  D32  is  (Zc  +1,  ZN_b  _c^)  where  <  ...  <  ZN_fa  are  the  ordered  elements 


of  D-2  and  C?  *  (a/2)  (N  -  bJA'“'/2  +  (H  -  b2  -  l)/2.  Now  let 


1/2  , 


32 


V  *  X  WC1+1  +  ZC2+1 


V  ^  A  WN-Cl-bl  +  (1  '  X)  ZN-C2-b2 


(4.2) 


where  0  <  X  <  1.  Then  we  propose  (9  *,  9*)  as  a  confidence  interval  for 
""  L  U 

*  ic  A  A 

Theorem  5.  (0^  ,  9y  )  is  asymptotically  a  (1  -  a)  100%  confidence  interval 

A  ^  A  ^ 

for  0;  ie,  lim  P(0T  <  0  <  0T  )  =  1  -  a. 

Ld  U 

N-wo 


0. 


Proof:  As  N  tends  to  infinity,  we  have  P(9  <  Wc  +^) ,  P(0  <  Z c  +^) , 

P(0  >  w*t  u  „  )  and  P(0  >  Z..  ,  _  )  each  tending  to  a/2.  But 

N-b^-C^  n-b2~C2 

min  (WCi+1,  ZC2+1)  <  0L*  <  max  {W^,  Z^}  and 
^  {WN-b1-C1*  ZN-b2-C2}  -9!-DaX  {WN-brC2‘  Vb2-C2}‘ 

A  A  ^  * 

Combining  these  results  implies  that  P(9  <  0L  )  tends  to  a/2  and  P(0  >  0y  ) 
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tends  to  a/2  which  proves  the  theorem. 

Next,  we  investigate  the  efficiency  of  this  interval  with 
1  *  (1  -  Y^)  /  (2  -  ”  Y2).  so  that  the  weights  are  proportional  to  the 

number  of  elements  in  the  sets  and  Notice  that,  in  this  case, 


9L  “  X  (XC1+1  +  '  2  +  (1  "  X>  V1  +  XC2^2+1)/2 


(4.3) 


9U  ”  X  ^-bj-Cj.  +  XN-C1)  /  2  +  d  ~  X>  (XN-b2-C2  +  XN-C2)  /  2 
A  A  * 

Thus,  (9y  -  9^  )  is  just  a  linear  combination  of  the  lengths  of  the  con¬ 

fidence  intervals  based  on  and  D^2,  respectively.  As  a  result,  we  can 
use  earlier  results  on  efficiencies  to  compute  the  efficiency  of  the  interval 

^  it  A  it 

(0L  ,  9y  ) .  In  particular,  it  can  be  shown  that 


~  *  *  * 

1/2  ®tj  ”  9L  2 

[n1/2  ( g  ■  ■■  - — )  r 

2  (y/2) 


(4.4) 


in  probability. 

Using  (4.4),  Table  6  compares  this  interval  with  other  commonly  used  confidence 

/V  A 

intervals  for  9  as  well  as  with  (0L>  8  )  with  Y^  *  *2,  Y2  *  *6*  We  see  that 
considerable  efficiency  is  lost  for  both  the  logistic  and  double  exponential 

A  ^  A  ^  A  A 

distributions  when  (9T  ,  9  )  replaces  (9  ,  9  ).  This  is  in  contrast  to 

Li  U  L  U 

the  estimation  case  where  the  efficiency  of  9  closely  paralleled  that  of 

/A 

9.  Thus,  we  see  that  in  the  case  of  interval  estimation,  a  somewhat  higher 
price  in  terms  of  efficiency  must  be  paid  to  get  the  additional  simplicity 

A  ^  A  ^  /\  A 

of  (9^  ,  0^  )  over  (9^,  9^) .  Once  again,  however,  both  of  these  methods 
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perform  quite  well  at  the  normal  model,  and  can  be  superior  to  the  interval 
based  on  the  sample  mean  for  non-normal  models. 

-  Table  6  about  here  - 
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5.  Extensions  and  Conclusions 


First,  we  briefly  summarize  without  proofs,  the  extensions  of  the 
previous  results  to  the  two-sample  location  problem.  In  particular,  we 
point  out  that  the  scores  function  as  a  step  function  produces  a  two-sample 
rank  statistic  and  corresponding  point  and  interval  estimators  with  the 
same  asymptotic  relative  efficiencies  as  those  considered  in  the  one-sample 
problem.  For  a  complete  discussion,  see  Markowski  (1980). 

Let  X.  <  ...  <  X  and  Y,  <  ...  <  Y  be  the  ordered  values  of  random 
i  mi  n 

samples  from  distributions  with  distribution  functions  F(x)  and  F(y  -  6), 

respectively.  Suppose  without  loss  of  generality  m  <  n. 

k  i  -  0 . [(1-Yk)  N/2] 

;  (5.1) 

0  i  -  [ (1~Y1)  N/2]  +  1 . [(1+Yl)  N/2) 

k  i  *  [ (1+Y^)  N/2]  +  1,  ...,  N 


Let  (^.(i)  =  j 


where  0  _<  y^  <  Y  2  <  •••  <  Yk  .£  1»  k  <  N  =*  m  +  n,  k  fixed. 

n 

Consider  SN  ^(6)  -  Z  c^  (1^(5))  where  11^(5)  is  the  rank  of  Y^  -  5  among 
’  i»l 

X^,  ...,  Xm,  Y^  -  6,  . ..,  Yn  -  6.  Now  with  Yg  “  0  and  Y  *  1  let 


8t  -  t (1  +  Yt)  N/2],  g_t  =  [(1  -  Yt)  N/2]  t  -  0 . k  +  1. 

k 


(5.2) 


Theorem  6.  If  5  *  0,  then  S  ,  (0) 


I  i  W.  where 
i—k 


P(W±  -  t) 


gi  “  *i-lx  ,N  "  (gi  “  gl-l> 


L)  ( 


n-t 


)/  O  for  t  -  0,  1, 


min  (n,  g±  -  g^L  In  addition,  E(SN>  k(0))  *  0,  Var  (S^  k(0)) 
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k  - 

2  mn  (  Z  i4  (g^+^  ~  g^))  /  N(N  -  1),  and  for  k  fixed,  the  limiting  distribu- 

1/2 

tion  (as  N  -*■  ®)  of  SN  k(0)  /  (Var  SN  k(0))  '  is  N(0,  1). 

Next,  we  consider  the  H-L  estimator  derived  from  ^(5),  which  will 

A 

be  denoted  by  6^. 

Theorem  7.  Let  gt  be  as  in  (5.2),  and  -  {Y.  -  :  i  +  j  *  g  +  1, 

/.  J 

t  ■  +  1,  ...  +  k}.  Then  6^  =  med  L^.  In  particular, 

A 

<5^  *  med  ■  med  (P^  U  P_^} 

(5.3) 


6 2  -  med  L2  -  med  {Pj.  U  P  U  P_2  U  P21 
where  P£  »  {Y^  -  Xi  »  i  +  j  *  gt  +  l). 

Notice  that  the  elements  within  any  set  P  are  automatically  ordered  and  so 

A  /v 

computation  of  5^  required  the  same  steps  as  computation  of  for  which 
an  algorithm  was  previously  given.  Computation  of  6^  requires  the  determina¬ 
tion  of  the  median  of  the  union  of  four  ordered  sets.  The  algorithm  considered 
here  can  be  extended  in  a  straightforward  way  to  allow  simple  computation  of 


In  the  same  way,  it  can  be  shown  that  a  (1  -  a)  100%  confidence  interval 

A  A  A 

for  6  has  the  form  (6^,  6y)  where  6^  =  (C  +  1)  'st  smallest  element  in 

A 

6y  ■  (N^  -  C)  'th  smallest  element  in  (5.4) 

N/  ■  the  number  of  elements  in  L.  and  C  satisfies  P(S  .(0)  <  C)  =  a/2. 

X  K.  g.Q  N»  K 

Thus,  we  see  that  computation  of  this  interval  requires  the  same  type  of 
algorithm  as  is  needed  to  compute  the  estimator.  Also,  approximate  confidence 
intervals  can  be  obtained  using  the  normal  approximation  given  in  Theorem  6. 

In  addition,  all  the  efficiency  comparisons  in  the  one-sample  case  remain 
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valid  when  comparing  these  Inference  methods  with  other  two-sample  methods. 
Also,  the  simplicity  of  the  methods  remains  In  the  two-sample  problem.  Thus, 
once  again  we  have  proposed  a  class  of  inference  methods  which  are  desirable 
if  efficiency,  simplicity,  and  robustness  are  considered. 

Lastly,  the  scores  function  considered  can  be  adapted  to  the  problem 
of  making  inferences  about  the  slope  parameter  in  a  simple  linear  regression 
model  when  the  form  of  the  error  distribution  is  unspecified.  Once  again, 
the  resulting  inference  procedures  will  retain  the  efficiencies  of  those  which 
we  have  discussed. 
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APPENDIX 


Lemma  1.  Let  ft  be  the  class  of  symmetric,  unimodal  densities,  then 


inf  e(0(Y1),  X)  =  4  (1-Y-i)2* 
fefl  1  31 


(A.  1) 


Proof : 


2  -1  ^+Y1 

4r[Fi(-^)] 

e(e(Yj).X)  -  - -  0“ 


(A.  2) 


We  assume  WLOG  that  f(0)  =1.  We  wish  to  find  f  £  £2  which  minimizes 
-1  ^+Y1  2 

f  [F  (-y^)]  for  fixed  Yr  Let 


fc(x) 


°  <  x  <  Yl/2 


(A.  3) 


—  Yx/2  <  x  <  c  +  Yj_/2 


and  f  (-x)  =  f  (x) . 
c  c 

Then  it  is  easily  verified  that  f  is  a  density  function  and  it  is  clear 

that  f  may  be  approximated  arbitrarily  closely  by  a  density  f^  in  £2. 

1+Y  1-Y,  j  ,  Yi  o  1-Y,  2 

Furthermore,  f c [Fc  (— j — )]  *  2c~  and  af  =  *3  ^  2~  ^  +  — 3 — 

c 

Yl(1-Yl)  Y12(1'Y1} 

- 2  c  +  4  • 


-1  ^+"^1  2 

Now  for  f[F  (  — - —  )]  fixed,  0^  is  minimized  for  unimodal  densities  by 


choosing  f  to  be  uniform  for  x  Y^/2.  Hence,  the  density  f  will  be 
made  least  favorable  by  letting  c  -*■  00 .  In  this  case. 


e(0(Y ,),X)-  -Ar- 
rl 


2c 


a 


2 

f 


and  letting  c  -*•  00  we  have 


inf  e(0(Y1),X)  - 
f  eQ 


_4 _ 

^1 


)(■ 


1-Y 

~F" 


l 

3 


U-Yi>2  • 


Finally,  we  note  that  f£  £  F^  but  since  it  can  be  approximated  arbitrarily 
closely  by  f.^  e  F^,  the  result  remains  valid. 

Notice  that  the  least  favorable  density  is  1  until  the  point  —  and 
then  becomes  arbitrarily  small  and  uniform  in  the  tails.  We  now  turn 

A 

to  the  result  comparing  0(Y^,Y2)  and  X.  First,  we  need  a  preliminary 
lemma. 

Lemma  2 . 

A  A 

lim  ef  (0(Y1,Y2),X)  <_  ef  (SCY^Yj)  .X) 
c-*<x>  c 


for  any  f  £  where  f£  is  given  by  A. 3. 

Proof:  Again  take  f(0)  =  1.  It  is  clear  from  the  proof  of  Lemma  1 


that  the  least  favorable  density  must  be  of  the  form: 


1 


0  <  x  <  Y1/2 


g(x) 


a  Yj^/2  £  x  <  y2/2 


Y2/2  <  x  <  Y2/2  +  d 


where  b  <_  a  <  1  and  d  _>  0. 

But  by  the  previous  lemma,  for  fixed  a,  e(6(YpY2)  >X)  minimized  by 
letting  b  -*•  0.  Furthermore,  for  any  fixed  a  and  b,  we  may  decrease  the 
efficiency  by  choosing  a'  <  a  and  letting  b  -*■  0.  Thus,  the  least  favorable 
choice  for  g  makes  a  arbitrarily  small  and  b  *  0.  But  this  is  just  the 
density  f  with  c  becoming  arbitrarily  large,  proving  the  result. 

Proof  of  Theorem  4. 

A  _  A  __ 

By  Lemma  2,  inf  e(9(Y, ,Yo) ,X)  *  lim  ef  (9(Y, ,Y,),X) 
fefl  c-*»  c 


1“Y  ^  1~Y  ^  2 

( - -  + - -)*  2 

-  lim  v  2c  2c  ‘  af 

^  1-Y2  +  i(Y2_Yl) 

(1~Y1) 2  (1-Y1)3 

“  i  •  3  ■  i 

!"y2  +  4^2"Y1J  3(1-Y2  +  ^(Y2-Y1)) 


Table  1  Maximum  Efficacy  of  T. 


Kurtosis 

c2* 

n 

Yi 

e(Yi, 

3 

.81 

.46 

.81 

4.2 

.296 

1/3 

o\ 

00 

7.65 

.74 

.38 

.84 

6 

1 

0 

1 

C  is  maximum  squared  efficacy  of  TM 

w  > 

*  * 
is  value  of  which  maximizes  C 


e(TH  LRT)  is  efficiency  of  T  .  with  y.  -  y,  to  likelihood  ratio 

N,  1|  A  ^  ^ 

test  (LRT) 


2.  Efficacy  of  2  as  a  function  of  Y^»  Y2  and 

c2 


1 

*2 

A3 

A4 

* 

637 

.250 

.595 

1 

,742 

.288 

.692 

.876 

,819 

.310 

.757 

.833 

.870 

.316 

.795 

.805 

.870 

.304 

.780 

.807 

.829 

.312 

.765 

.865 

.888 

* 

.320 

.808* 

.720 

.902 

.311 

.807 

.714 

.876 

.304 

.789 

.556 

* 

.912 

.300 

.803 

.533 

.903 

.295 

.764 

.544 

.866 

.283 

.775 

.457 

.884 

.280 

.774 

.445 

.885 

.276 

.746 

.450 

.830 

.250 

.713 

.360 

.688 

.357 

[I, 


3.  Efficiency  of  TN  2  with  Y^  ■  .2,  y2  ■  .6 


A2 

A3 

\ 

'<T», 

2,  S) 

1.40 

1.28 

1.32 

» 

•<V 

2,  W) 

.930 

.96 

.99 

• 

e<Y 

2,  ^ 

.890 

1.05 

1.13 

1. 

2,  T  N, 

.976 

1 

1 

• 

2,  LRT) 

.890 

.961 

.960 

• 

S  -  sign  test,  W  -  Wilcoxon,  t  -  t-test,  TN  2  "  optimal  choice 


Table  4 


ARE  (9 
ARE  (9 
ARE  (6 
ARE(0* 

Note: 


ARE  of 

CD  > 

» 

• 

2,  .6) 

with  respect 

to  other  estimators 

o 

t-h 

a> 

A1 

*2 

A3 

A4 

X) 

.87 

1.04 

1.12 

1.44 

HL) 

.91 

.96 

.98 

.96 

X) 

1.37 

1.28 

1.30 

.72 

0(.2, 

.6)) 

.98 

.99 

.99 

1.00 

X  -  sample  mean,  HL  -  median  of  Walsh  averages,  X  -  sample  median 


A  * 

6.  ARE  of  (0L  , 

^  ft 

0y  )  relative 
(Yi-  ■ 

to  other 

.2,  Y2  -  • 

interval  estimators 

6) 

of  ( 

A1 

A2 

A3 

A4 

ASE(1,  2) 

.87 

.91 

1.10 

1.24 

ARE(1,  3) 

1.37 

1.11 

1.28 

.62 

AKE(1,  4) 

.91 

.84 

.96 

.83 

ARE(1,  5) 

.98 

.87 

.97 

.86 

A  ^  If 

1  -  (0L  (.2,  .6),  Qy  (.2,  .6)),  2  -  t-interval,  3  -  sign  test  interval 

A  A 

4  -  Wilcoxon  test  interval,  5  -  (0y(.2,  .6),  9y(.2,  .6)) 
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